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ABSTRACT 



Context. One of the typical features shown by observations of solar prominence oscillations is that they are damped in time and 
that the values of the damping times are usually between one and three times the corresponding oscillatory period. However, the 
mechanism responsible for the attenuation is still not well-known. 
Aims. Thermal conduction, optically thin or thick radiation and heating are taken into account in the energy equation, and their role 
on the attenuation of prominence oscillations is evaluated. 

Methods. The dispersion relation for linear non-adiabatic magnetoacoustic waves is derived considering an equilibrium made of 
a prominence plasma slab embedded in an unbounded corona. The magnetic field is orientated along the direction parallel to the 
slab axis and has the same strength in all regions. By solving the dispersion relation for a fixed wavenumber, a complex oscillatory 
C/3 i frequency is obtained, and the period and the damping time are computed. 

Results. The effect of conduction and radiation losses is different for each magnetoacoustic mode and depends on the wavenumber. In 
the observed range of wavelengths the internal slow mode is attenuated by radiation from the prominence plasma, the fast mode by the 
combination of prominence radiation and coronal conduction and the external slow mode by coronal conduction. The consideration of 
the external corona is of paramount importance in the case of the fast and external slow modes, whereas it does not affect the internal 
slow modes at all. When a thinner slab representing a filament thread is considered the fast mode is less attenuatted whereas both 
internal and external slow modes are not affected. 

Conclusions. Non-adiabatic effects are efficient damping mechanisms for magnetoacoustic modes, and the values of the obtained 
damping times are compatible with those observed. 
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O 1. Introduction oscillations (P > 40 min), although very short-periods of less 

kJ than 1 min (Balthasar et al. 1993) and extreme ultra-long-periods 

Prominences are dense coronal structures which appear as thin, of more man 8 hours (Fou i lon et aL 2 004) have been reported. 
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dark filaments on the solar disc when observed in Ha. On the Nevertheless, the value of the period seems not to be related with 

^ ■ contrary, they show up as bright objects above the solar limb. the nature or the sourc£ of the and probably is linked to 

Ctf ; The coronal magnetic field is responsible for the support of the promin ence eigenmode that is excited. There are also a few 

prominences against gravity, and it also plays a fundamental role determinations of the wavelength and phase speed of standing 

in the thermal confinement of the cool prominence plasma em- oscillations and pr0 pagating waves in large regions of promi- 

bedded in the much hotter coronal environment. Nevertheless, nences ( M olowny-Horas et al. 1997; Terradas et al. 2002) and 

the structure, orientation and strength of the magnetic field in in single filamen t threads (Lin et al. 2007). On the other hand, 

prominences and the surrounding corona is still enigmatic and several obser vations (Molowny-Horas et al. 1999; Terradas et al. 

not well-known. High resolution observations reveal that promt- 2002) haye informed abo ut the evidence of the attenuation of the 

nences are composed by numerous very thin, thread-like struc- oscillations in Doppler velocity time series, which is a common 

tures, called fibrils, piled up to form the body of the prominence feature obserV ed in large areas. By fitting a sinusoidal function 

(Lin et al. 2003; Lin et al. 2005, Lin et al. 2007) and measures multiplie d by a factor exp(-f/T D ) to the Doppler series, these 

also indicate that magnetic field lines are orientated along these authors have obtaine d values of the damping time, t d , which 

thin threads. are usua jjy b etween 1 and 3 times the corresponding oscilla- 

The observational evidence of small-amplitude oscillations tory period. The reader is referred to some recent reviews for 

in quiescent solar prominences goes back to 40 years ago more information about the observational background (Oliver & 

(Harvey 1969). The amplitude of these oscillations typically Ballester 2002, Wiehr 2004, Engvold 2004, Ballester 2006). 
goes from less than 0. 1 km s 1 to 2-3 km s 1 , and have been 

historically classified, according to their periods, in short- (P < From the theoretical point of view, small-amplitude promi- 
10 min), intermediate- (10 min < P < 40 min) and long-period nence oscillations can be interpreted in terms of linear magne- 
tohydrodynamic (MHD) waves. Although, there is a wide bib- 



Send offprint requests to: R. Soler liography of works that investigate the ideal MHD wave modes 



2 



R. Soler et al.: The effect of the solar corona on the attenuation of prominence oscillations I 



supported by prominence models (see Oliver & Ballester 2002 
for an extensive review of theoretical studies), the investigation 
of the wave damping has been broached in few papers. By re- 
moving the ideal assumption and including dissipative terms in 
the basic MHD equations, several works have studied the atten- 
uation of prominence oscillations considering radiative losses 
based on the Newtonian law of cooling with a constant relax- 
ation time (Terradas et al 2001), or performing a more complete 
treatment of non-adiabatic effects, assuming optically thin radi- 
ation, heating and thermal conduction (Carbonell et al. 2004; 
Terradas et al. 2005). The main conclusion of these previous 
studies is that only the slow wave is attenuated by thermal 
effects, radiation being the dominant damping mechanism in 
the range of typically observed wavelengths in prominences, 
but the fast wave remains practically undamped. On the other 
hand, Forteza et al. (2007) proposed ion-neutral collisions as a 
damping mechanism on the basis that prominences are partially 
ionised plasmas, but they found that this mechanism is only ef- 
ficient in attenuating the fast mode in quasi-neutral plasmas, the 
slow mode being almost unaffected. 

In the light of these referred studies, it is likely that non- 
adiabatic effects are the best candidates for the damping of 
small-amplitude oscillations, at least for slow modes. However, 
previous results do not asses the influence of the corona. The 
main aim of the present work is to perform a step forward in the 
investigation of the effect of non-adiabatic mechanisms (radia- 
tion losses, thermal conduction and heating) on the time damp- 
ing of prominence oscillations. We consider a slab model with 
a longitudinal magnetic field and take into account the external 
coronal medium. So, we explore for the first time the joint ef- 
fect of prominence and coronal mechanisms on the attenuation 
of oscillations. The magnetoacoustic normal modes of this equi- 
librium have been previously investigated by Edwin & Roberts 
(1982) and Joarder & Roberts (1992) in the adiabatic case. Later, 
a revision of these works has been done in Soler et al. (2007), 
hereafter Paper I, and the normal modes have been studied and 
reclassified according to their magnetoacoustic properties. 

This paper is organised as follows. The description of the 
equilibrium model and the linear non-adiabatic wave equations 
are given in Sect. [2] whereas the dispersion relation for the mag- 
netoacoustic modes is derived in Sect. [3] Then, the results are 
plotted and investigated in Sect. [4] Finally, Sect. [5] contains the 
conclusions of this work. 



2. Equilibrium and basic equations 

Our equilibrium configuration (Fig. [TJ is made of a homoge- 
neous plasma layer with prominence conditions (density p p and 
temperature T p ) embedded in an unbounded corona (density p c 
and temperature T c ). The coronal density is computed by fix- 
ing the coronal temperature and imposing pressure continuity 
across the interfaces. The magnetic field is Bo = Boe x , with Bo 
a constant everywhere. Both media are unlimited in the x- and 
y-directions. The half-width of the prominence slab is z p . 

The basic magnetohydrodynamic equations for the discus- 
sion of non-adiabatic processes are: 



Dp 

Dt H 

Dv 1 
Dt p 



(1) 



(2) 



2l-IP£E +(y -l )[pL Q 3> T)-V-(K-VT)]=0, (3) 
Dt p Dt 



A 



Corona 



Prominence 



2z„ 



Corona 



Bo > T € , p c 



Fig. 1. Sketch of the equilibrium. 

Table 1. Parameter values of the radiative loss function 
corresponding to the considered regimes. The three promi- 
nence regimes represent different plasma optical thicknesses. 
Prominence (1) regime corresponds to an optically thin plasma 
whereas Prominence (2) and Prominence (3) regimes represent 
greater optical thicknesses. All quantities are expressed in MKS 
units. 



Regime 


x" 


a 


Reference 


Prominence (1) 


1.76 x irr lj 


7.4 


Hildner(1974) 


Prominence (2) 


1.76 x 1CT 53 


17.4 


Milne et al. (1979) 


Prominence (3) 


7.01 x lfr 104 


30 


Rosneretal. (1978) 


Corona 


1.97 x 10 24 


-1 


Hildner(1974) 



y=Vx(vx8), 
at 

V ■ B = 0, 

pRT 
P = —=-, 
P 



(4) 
(5) 
(6) 



where tt = 4 + v 

Dt at 



V is the material derivative for time vari- 
ations following the motion and all quantities have their usual 
meaning. Equation (f3]i is the energy equation, which in the 
present form takes into account non-adiabatic effects (radiation 
losses, thermal conduction and heating) and whose terms are ex- 
plained in detail in Carbonell et al. (2004) and Terradas et al. 
(2005). Following these works, only thermal conduction par- 
allel to the magnetic field is assumed and we use the typical 
value for the parallel conductivity in prominence and coronal 
applications, k\\ = 10 _11 r 5 ^ 2 WnT 1 K _1 . Radiative losses and 
heating are evaluated together through the heat-loss function, 
L(p, T) = x*pT a - hp a T h , where radiation is parametrised with 
X* and a (see Table [TJ and the heating scenario is given by ex- 
ponents a and b. The heating mechanisms taken into account in 
this work are (Rosner et al. |1978t Dahlburg & Mariska [T988] ): 

- constant heating per unit volume {a - b - 0); 

- constant heating per unit mass (a — I, b = 0); 

- heating by coronal current dissipation (a — b — 1); 

- heating by Alfven mode/mode conversion (a-b- 7/6); 

- heating by Alfven mode/anomalous conduction damping 
(a= l/2,b= -1/2). 

Following the same process as in Carbonell et al. (2004), 
we consider small perturbations from the equilibrium state, lin- 
earise the basic Eqs. (UJ)-© and obtain their Eqs. (9)-(14). Since 
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our model is unlimited in the x- and y-directions, we assume all 
perturbations are in the form f\ (z) exp i((ot + k x x + k y y), and con- 
sidering only motions and propagation in the xz-plane (v v = 0, 
k y = 0), which excludes Alfven waves, the linearised equations 



become 



icopi + po Uk x v x + -^1 = 0, 



iu>pov x - -ik x p\, 
dpi 



ILOPQV, = - 



Bo ., dB lx 

+ — \ik K B\, 

dz p \ dz 



(7) 
(8) 
(9) 



icoipi -c 2 pi) = -(y- l)(^/sr l |7'i + —to p p\ + ^-cottA, (10) 
x \ Po la I 



dz' 

i<jjB\ z — Boik x v z , 

where c 2 — — is the adiabatic sound speed squared and 



(11) 
(12) 



co p = — (L +pqL„) , co T = —ToL T , 
po x ' po 

L p , Lj being the partial derivatives of the heat-loss function with 
respect to density and temperature, respectively, 



dT 



Now, it is possible to eliminate all perturbations in favour of 
v z to obtain a single differential equation 



d 2 v. 



dz 2 
in which 

(cJ-ky^-kiK 1 ) 



K = 



(v 2 A + A 2 )(co 2 -k 2 x c 2 T ) 



(13) 



(14) 



where v\ 



^ is the Alfven speed squared. A 2 and c 2 , are the 
modified sound and cusp (or tube) speed squared, respectively, 

(y - 1) (^11^ w p) + iju 



A z = -± 

r 



,2- W 

T ~v 2 +A 2 ' 



(15) 



(16) 



Expressions for the perturbations in terms of v z are given in 
App. A. In all the following formulae, subscripts p or c denote 
quantities computed using prominence or coronal values, respec- 
tively. 



3. Dispersion relation 

We impose some restrictions on the solutions of Eq. (13[ in or- 
der to obtain the dispersion relation for the linear non-adiabatic 
magnetoacoustic waves. We restrict this analysis to body waves 
which are evanescent in the corona, since we are looking for 
solutions which are essentially confined to the slab. For such so- 
lutions, v z (z) is of the form 



v z (z) = 



Aiexp^c^ + Zp)], 

A 2 cos (kzpz) + A 3 sin (k^z) , 

A 4 exp [-k zc (z - Zp)J , 



if z < -z p , 

if -z p < z < z p , 

if Z > Zn. 



(17) 



with %(k zp ) > and %(k zc ) > 0. 

Imposing continuity of v z and the total (gas plus magnetic) 
pressure perturbation across the interfaces, we find four alge- 
braic relations between the constants A\, A 2 , A3 and A4. The 
non-trivial solution of this system gives us the dispersion rela- 
tion 

g {ky Ac - co 2 ) k zp { 2 } M ± {k 2 A P - - 2 ) h c =0, (18) 



where cot/tan terms and + signs are related with the symmetry of 
the perturbations. The cot term and the + sign correspond to kink 
modes (A3 = 0), whereas the tan term and the - sign correspond 
to sausage modes (A2 = 0). 

The dispersion relation for the magnetoacoustic waves pre- 
sented in Eq. ( fT8l is equivalent to the relation investigated in 
Edwin & Roberts (1982) and Joarder & Roberts (1992), and re- 
vised Paper I, in the case of adiabatic perturbations, since all 
non-adiabatic terms are now enclosed in k zp and k, c through 
Eq. 01- 



4. Results 

Now, we assume Prominence (1) conditions inside the slab (i.e. 
an optically thin prominence) and a heating mechanism given 
by a = b = 0. Unless otherwise stated, the following equilibrium 
parameters are considered in all computations: T p = 8000 K, 
p p = 5 x 10- 11 kg trr 3 , T c = 10 6 K, p c = 2.5 x 10~ 13 kg trr 3 , 
Bo = 5 G andzp = 3000 km. The solution of the dispersion 
relation (Eq. lfT8l ) for a fixed real k x gives us a complex fre- 
quency co — wr + ico\. We then compute the oscillatory period, 
the damping time and the ratio of the damping time to the pe- 
riod because this is an important quantity from the observational 
point of view, 



2n 

P= — , 

COr 



1 

td = — , 



td 
P 



1 

2n a>i 



Since we are interested in studying the behaviour of the most 
relevant solutions of the dispersion relation, we only compute the 
results for the fundamental modes, which are labelled, accord- 
ing to the classification of Paper I, as internal and external slow 
modes and fast modes. The band structure described in Paper I is 
slightly modified when non-adiabatic terms are considered (see 
Fig. |21). The phase speed of the internal slow modes is now en- 
closed in the band %{cj p ) < co R /k x < %(A p ). The adiabatic fast 
modes exist in two separated bands in the phase speed diagram 
due to the presence of a forbidden region {cj c < coR/k x < c sc ), 
but now the forbidden band is avoided and a continuous fast 
mode is found with va p < tOR/k x < va c - Finally, and like in 
the adiabatic case, among the external slow modes only the fun- 
damental kink one exists as a non-leaky solution in a restricted 
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wavenumber range and couples with the fundamental fast kink 
mode. Its phase speed is u>n/k x « !R(A C ). Therefore, we see that 
in the non-adiabatic case A plays the role of c s in the adiabatic 
case. 




Fig. 2. Phase speed versus the dimensionless wavenumber for 
the three fundamental oscillatory modes. Solid lines denote kink 
modes whereas dotted lines correspond to sausage modes. The 
shaded zones are projections of the forbidden (or leaky) regions 
on the plane of this diagram. Note that the vertical axis is not 
drawn to scale. 

In Fig.[3]P, tq and Tu/P are represented for the fundamental 
modes and for a range of the longitudinal wavenumber between 
lO- 10 !^ 1 and 10 3 m" 1 . The shaded zones correspond to wave- 
lengths between 5 x 10 3 km and 10 5 km, the typically observed 
values. It turns out that the values of the period are very sim- 
ilar to those obtained in the adiabatic case (Joarder & Roberts 
1992, Paper I). The damping time presents a strong dependence 
with the wavenumber and its behaviour is very different from 
one mode to another. This fact suggests that the non-adiabatic 
mechanisms can affect each mode in a different way (Carbonell 
et al. 2004). This is studied in detail in Sect 14.11 Observations 
show that prominence oscillations are typically attenuated in a 
few periods (Terradas et al. 2002), so a damping time of the order 
of the period is expected. In our results, the fundamental modes 
present values of td /P in the range 1 to 10 in the observed wave- 
length region, which is in agreement with observations. 

4. 1 . Regions of dominance of the damping mechanisms 

The importance of the different non-adiabatic terms included in 
the energy equation (Eq. ]3)) depends on the wavenumber. In 
order to know which is the range of dominance of each mech- 
anism, we compare the damping time obtained when consider- 
ing all non-adiabatic terms (displayed in the middle column of 
Fig. [3j with the results obtained when a specific mechanism is 
removed from the energy equation. With this analysis, we are 
able to know where the omitted mechanism has an appreciable 
effect on the damping. The results of these computations for the 
fundamental kink modes (Fig. |4j are summarised as follows: 

- The fundamental internal slow kink mode is not affected by 
the mechanisms related with the corona. This is a conse- 
quence of the nature of this mode, which propagates strictly 
along the prominence without disturbing the corona (see Fig. 



4, top row, of Paper I). For this reason, in the adiabatic case 
it is also independent of the coronal conditions. On the other 
hand, the prominence-related mechanisms show different ef- 
fects in two different ranges of k x . For k x < 10~ 3 trT 1 promi- 
nence radiation dominates, while for k x > 10~ 3 m _I promi- 
nence conduction is the dominant mechanism. Beginning 
from small values of the wavenumber, prominence radiation 
becomes more efficient as k x grows and the damping time 
falls following a power law until k x « 10 m _1 , where td 
saturates in a plateau between k x « 10~ 5 and k x « 10~ 3 trT 1 . 
Then, prominence conduction becomes the dominant mech- 
anism and the damping time falls again until k x * 10 _1 rrT 1 
where a new plateau begins. This last part of the curve cor- 
responds to the isothermal or superconductive regime, in 
which the amplitude of the temperature perturbation drops 
dramatically (Carbonell et al. 2006). Prominence radiation 
is responsible for the attenuation of the slow mode in the 
observed wavelength range. An approximate dispersion re- 
lation for the internal slow modes is included in App. [B] 

- The fundamental fast kink mode is affected by the four 
mechanisms. For k x < 3 x 10~ 9 trr' coronal radiation dom- 
inates but for 3 x lO^m -1 < k x < 5 x lO^m -1 the effect 
of coronal conduction grows and becomes the main damping 
mechanism. Then, for k x > 5xlO _7 m _I the corona loses dra- 
matically its influence and prominence mechanisms become 
responsible for the attenuation of this mode. First, promi- 
nence radiation is dominant in the range 5 x 10~ 7 m _I < 
k x < 10 _3 m , then prominence conduction governs the 
wave damping for k x > 10 _3 m _1 and finally the isother- 
mal regime begins for k x « 10° m -1 . The minimum of Td 
occurs into the coronal conduction regime, for the value of 
k x which corresponds to the coupling with the external slow 
mode. The transition between the coronal conduction regime 
and the prominence radiation regime occurs in the observed 
wavelength range. The reason for the sensitivity of the fast 
mode damping time on prominence and coronal conditions 
is that this wave has a considerable amplitude both inside the 
prominence and in the corona, the later becoming more im- 
portant for long wavelengths (see the second and third rows 
of Fig. 4 of Paper I). 

- The behaviour of the damping time of the fundamental exter- 
nal slow kink mode is entirely dominated by coronal mech- 
anisms whereas the prominence mechanisms do not affect it 
at all. This behaviour is a result of the negligible amplitude 
of this wave in the prominence (see the fourth and fifth rows 
of Fig. 4 of Paper I). For k x < 3 x 10~ 9 trT 1 coronal radiation 
dominates, but for shorter wavelengths coronal conduction 
becomes more relevant and is responsible for the damping 
in the observed wavelength range until the frequency cut-off 
is reached. At the cut-off, tq has a value of the order of the 
period. 



Regarding the fundamental sausage modes, the behaviour of 
the internal slow sausage mode is exactly that of the slow kink 
mode, so no additional comments are needed. The fundamen- 
tal fast sausage mode (Fig. ID presents the same scheme as the 
fundamental fast kink mode for k x > 10~ 8 m . The main dif- 
ference between the fast kink and sausage modes happens in the 
observed wavelength range, where the effect of coronal conduc- 
tion on the sausage mode is less relevant. If coronal conduction 
is omitted, the fundamental fast sausage mode is not able to tra- 
verse the forbidden region in the dispersion diagram and then 
shows frequency cut-offs as in the adiabatic case. This means 
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Fig. 3. Period (left), damping time (centre) and ratio of the damping time to the period (right) versus the longitudinal wavenumber 
for the fundamental oscillatory modes. Upper panels: internal slow kink (solid line), fast kink (dotted line) and external slow kink 
(dashed line). Lower panels: internal slow sausage (solid line) and fast sausage (dotted line). Shaded zones correspond to those 
wavelengths typically observed. Note the cut-offs of the external slow kink mode and the fast sausage mode. Prominence (1) 
radiation conditions have been taken for the prominence plasma and the heating scenario is given by a = b = 0. 
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Fig. 4. Damping time versus the longitudinal wavenumber for the three fundamental kink oscillatory modes: internal slow (left), 
fast (centre) and external slow (right). Different linestyles represent the omitted mechanism: all mechanisms considered (solid line), 
prominence conduction eliminated (dotted line), prominence radiation eliminated (dashed line), coronal conduction eliminated (dot- 
dashed line) and coronal radiation eliminated (three dot-dashed line). Prominence (1) radiation conditions have been taken for the 
prominence plasma and the heating scenario is given by a = b = 0. 



that coronal conduction causes the fast mode to cross the forbid- 
den region in the dispersion diagram in the non-adiabatic case. 



(r-i)TsP^^ Ts 



ypo 



Trad 



(20) 



Approximate values of k x for which the transitions between where Ts is the sound ttavel time and ^ an(J ^ ^ character . 

regimes take place can be computed by following a process simi- istic conductive and radiative time scales. Taking t s = 27i/k*c s , 

lar to that in Carbonell et al. (2006). The thermal ratio, d, and the the yalue of k * for which the condition d = r is satisfi ed is 
radiation ratio, r, quantify the importance of thermal conduction 
and radiation, respectively (De Moortel & Hood 2004), 



k* = 2npo ■ 



a- 1 



«\\ 



(21) 



d = 



(y - i)K\\T a p Q l r s 



7 2 PIT S 1 Tcond ' 



(19) 



Now, we use prominence values to compute k* for the promi- 
radiation-prominence conduction transition (fc*), and 



nence 
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Fig. 6. Damping time versus the longitudinal wavenumber for the fundamental internal slow kink mode (left), the fundamental fast 
kink mode (centre) and the fundamental external slow kink mode (right). The different linestyles represent different values of the 
prominence temperature: T p = 8000 K (solid line), T p = 5000 K (dotted line) and T p = 13000 K (dashed line). The heating scenario 
is given by a = b — and the optical thickness for the prominence plasma is Prominence (1). 
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Fig. 8. Same as Fig.|6]with B = 5 G (solid line), B = 2 G (dotted line) and B Q = 10 G (dashed line). 



coronal values for the coronal radiation-coronal conduction 



transition (fc*). This gives the values K * 1,7x10 m , and 
k* « 2.2 x 10~ 8 m . For the transition of the fast kink mode 
between the coronal conduction and the prominence radiation 
regimes, the boundary wavenumber k*^ c can be roughly calcu- 
lated by imposing d c = r p , that gives 



27rp p 



(22) 



and whose numerical value is k*^ 
wavenumbers for the transitions between different regimes are 
independent of the wave type, be it fast or slow, internal or ex- 



1.4 X 10~ 6 m" 1 . All these 



ternal (this agrees with Figs.|4]and|5j. On the other hand, the be- 
ginning of the isothermal regime can be estimated by following 
Porter et al. (1994). Considering Cj p /v^ p <k 1 and the approxi- 
mations wr * k x c sp for the slow wave and wr w & a va p for the 
fast wave, the critical wavenumber is 



2p p k B c 



^crit-slo\ 

K\\ p m p COS ft 

for the internal slow mode, and 
2p p k B v Ap 



^crit-fast 



K\\ p m p COS 2 fl' 



(23) 



(24) 
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Fig. 5. Same as Fig.|4]for the fundamental fast sausage mode. 

for the fast mode, where m p is the proton mass, ks is the 
Boltzmann constant and 9 is the angle between B and k. Taking 
cosi9 = 1 for simplicity, the approximate critical values are 
kcrit-siow ~ 1.7 x KT 1 nr 1 and ^ C1 -it-f ast * 9.1 x 1CT 1 m _1 . We 
note that all these approximate values describe correctly the tran- 
sitions between the diverse regimes shown in Figs. @] and |3J but 
their numerical values overestimate by almost an order of mag- 
nitude the actual critical wavenumbers. 

4.2. Exploring the parameter space 

4.2.1. Dependence on the equilibrium physical conditions 

In this section, we compute the solutions for different values 
of the equilibrium physical conditions. We only present the re- 
sults for the fundamental kink modes since they are equivalent to 
those of sausage modes. Figures|6]|7]and|8]display the damping 
time as function of k x for some selected values of the prominence 
temperature, the prominence density and the magnetic field, re- 
spectively. 

For the internal slow mode, a decrease of the prominence 
temperature or the prominence density raises the position of 
the radiative plateau and increases its length. The opposite be- 
haviour is seen when the density or the temperature are in- 
creased. However, the value of the magnetic field does not in- 
fluence the attenuation of this mode, such as expected for a slow 
wave. 

Increasing the value of the prominence temperature causes 
a vertical displacement of td of the fast mode in those regions 
in which prominence mechanisms dominate. The value of the 
prominence density has a smaller effect and its main influence 
is in changing the coupling point with the external slow mode, 
which moves to higher k x for greater values of the density. The 
magnetic field strength has a more complex effect on td and also 
modifies the coupling point. 

Finally, the external slow mode is only slightly affected by 
a modification of the prominence physical parameters since it 
is mainly dominated by coronal conditions, and the influence of 
the magnetic field is very small due to the slow-like magnetoa- 
coustic character of this solution. 

4.2.2. Dependence on the prominence optical thickness 

The optically thin radiation assumption is a reasonable approxi- 
mation in a plasma with coronal conditions but prominence plas- 
mas often are optically thick. In this section we compare the 
results obtained considering different optical thicknesses for the 



prominence plasma (see Fig.|9]for the fundamental kink modes). 
The results corresponding to the slow sausage mode have not 
been plotted since they are equivalent to those obtained for slow 
kink mode; those for the fundamental fast sausage mode, how- 
ever, are displayed in Fig.fTUl 

The variation of the prominence optical thickness modi- 
fies the prominence conduction-prominence radiation critical 
wavenumber, k* (see analytical approximation of Eq. EH ). For 
the internal slow mode, an increase in the optical thickness 
raises the position of the radiative plateau and shifts it to smaller 
wavenumbers. This fact causes an a priori surprising result in the 
observed wavelength range, since td has a smaller value for opti- 
cally thick radiation, Prominence (3), than for optically thin radi- 
ation, Prominence (1). Regarding fast modes, the damping time 
increases when the optical thickness is increased, but only in 
the region in which prominence radiation dominates. The value 
of td inside the observed wavelength range is partially affected 
and raises an order of magnitude for Prominence (3) conditions 
in comparison with the results for Prominence (1) conditions. 
Finally, the damping time of the external slow mode is not af- 
fected by the prominence optical thickness since it is entirely 
dominated by the corona, as it has been noticed in Sect. 14. II 




/<r s io~ s io~ 4 w~ 2 io° io s 

K (rn 1 ) 



Fig. 10. Same as Fig.|9]for the fundamental fast sausage mode. 



4.2.3. Dependence on the heating scenario 

Now, we compute the damping time for the five possible heating 
scenarios. For simplicity, we only consider the fundamental kink 
modes (Fig.fTTTi. Carbonell et al. (2004) showed that in a plasma 
with prominence conditions the different heating scenarios have 
no significant influence on the damping time. Nevertheless, in 
coronal conditions wave instabilities can appear depending on 
the heating mechanism. In our results, we see that the heating 
scenario affects the value of td only in the ranges of k x in which 
radiation is the dominant damping mechanism. The heating sce- 
nario has a negligible effect when prominence radiation domi- 
nates, since td is only slightly modified. On the contrary, wave 
instabilities appear in those regions in which coronal radiation 
dominates. Thermal destabilisation occurs when the imaginary 
part of the frequency becomes negative, so oscillations are not 
attenuated but amplified in time. Instabilities only occur in the 
fundamental fast kink and the external slow modes for very small 
values of k x , outside the observed wavelength range. 
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Fig. 9. Same as Fig.[6]with the prominence optical thickness given by Prominence (1) (solid line), Prominence (2) (dotted line) and 
Prominence (3) (dashed line) conditions. 
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Fig. 11. Same as Fig.|6]with the heating scenario given by a — b — (solid line); a = 1, b = (dotted line); a = b = 1 (dashed line); 
a — b — 7/6 (dot-dashed line); a = 1/2, b = — 1/2 (three dot-dashed line). 



4.3. Comparison with the solution for an isolated slab 

In order to assess the effects arising from the presence of two 
different media in the equilibrium, a comparison between the 
previous results and those corresponding to a single medium is 
suitable. So, we consider a simpler equilibrium made of an iso- 
lated prominence slab with the magnetic field parallel to its axis. 
The external medium is not taken into account. Magnetoacoustic 
non-adiabatic perturbations are governed by Eq. (fT3l l. and rigid 
boundary conditions for v z are imposed at the edges of the promi- 
nence slab, 

v z (-z p ) = v z ( Zp ) = 0. (25) 
Then, the solution is of the form 

v z (z) = Ci cos (kzpZ) + Ci sin (k zp zj , (26) 

and after imposing boundary conditions (Eq. ||25l ), we deduce 
the dispersion relation for the magnetoacoustic slow and fast 
non-adiabatic waves, 

kzpZ P = \n+^7i, (« = 0,1,2,...), (27) 
for the kink modes, and 

k zp z v = nn, (n = 1,2,3,...), (28) 

for the sausage modes. Inserting expressions (fT4l > and ( TT3T > for 
k~p and A p respectively, one can rewrite the dispersion relations 
(l27b and ( I28l i as polynomial equations in to. See App.lClfor the 
details. 



Next, considering only the fundamental kink modes for sim- 
plicity, we compute the period and the damping time and com- 
pare with those obtained when the surrounding corona is taken 
into account (Fig. fT2l . We see that there is a perfect agree- 
ment between both results in the case of the internal slow mode, 
whereas the solutions for the fast mode only coincide for inter- 
mediate and large wavenumbers, and show an absolutely differ- 
ent behaviour in the observed wavelength range and for smaller 
wavenumbers. Additionally, one must bear in mind that the ex- 
ternal slow mode exists because of the presence of the coronal 
medium, hence it is not supported by an isolated slab. 

In Paper I we proved that the internal slow mode is essen- 
tially confined within the prominence slab and that the effect of 
the corona on its oscillatory period can be neglected. Now, we 
see that the corona has no influence on the damping time either. 
On the other hand, the confinement of the fast mode is poor for 
small wavenumbers, the isolated slab approximation not being 
valid. As it has been noted in Section l4Tl the corona has an es- 
sential effect on the attenuation of the fast mode in the observed 
wavelength range. 

4.4. Application to a prominence fibril 

Since magnetic field lines are orientated along fibrils, our model 
can also be applied to study the oscillatory modes supported by 
a single prominence fibril. In order to perform this investigation, 
we reduce the slab half-width, z p , to a value according to the 
typical observed size of filament threads, which is between 0.2 
to 0.6 arcsec (Lin et al. 2005). Since these values are close to the 
resolution limit of present-day telescopes, it is likely that thinner 
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Fig. 12. Comparison between the solutions for a prominence 
plus corona system and for an isolated slab with prominence 
conditions. The upper panels correspond to the fundamental in- 
ternal slow kink mode and the lower panels to the fundamental 
fast kink mode. The solid lines are the solutions for a prominence 
plus corona equilibrium whereas the dotted lines with diamonds 
represent the solutions for an isolated slab. Prominence (1) pa- 
rameters and a = b = have been used in the computations. 



threads could exist. So, assuming now z p = 30 km, we compute 
P, To and T£,/P for the fundamental kink modes and compare 
these results with those obtained for z p = 3000 km. 

Such as displayed in Fig.[13j both internal and external slow 
modes are not affected by the width of the prominence slab since 
they are essentially polarised along the x-direction and so they 
are not influenced by the equilibrium structure in the z-direction. 
Nevertheless, the location of the cut-off of the external slow 
mode and the coupling point with the fast mode are shifted to 
larger values of k x when the slab width is reduced. On the other 
hand, the fast mode, which is responsible for transverse motions, 
is highly influenced by the value of z p . The Td curve for the 
fast mode is displaced to larger values of k x when smaller z p is 
considered. This causes that higher values of td/P are obtained 
in the observed wavelength range. Hence, these results suggest 
that local prominence oscillations related with transverse fast 
modes supported by a single fibril could be less affected by non- 
adiabatic mechanisms than global fast modes supported by the 
whole or large regions of the prominence. However, according to 
the results pointed out by Diaz et al. (2005) and Diaz & Roberts 
(2006), large groups of fibrils tend to oscillate together since the 
separation between individual fibrils is of the order of their thick- 
ness. In a very rough approximation one can consider that a thick 
prominence slab could represent many near threads which oscil- 
late together and that the larger the slab width, the more threads 
fit inside it. So, our results show that the slab size (i.e. the number 
of threads which oscillate together in this rough approximation) 
has important repercussions on the damping time of collective 
transverse oscillations, hence the oscillations could be more at- 
tenuated when the number of oscillating threads is larger. This 
affirmation should be verified by investigating the damping in 
multifibril models. 



5. Conclusions 

In this paper, we have studied the time damping of magne- 
toacoustic waves in a prominence-corona system considering 
non-adiabatic terms (thermal conduction, radiation losses and 
heating) in the energy equation. Small amplitude perturbations 
have been assumed, so the linearised non-adiabatic MHD equa- 
tions have been considered and the dispersion relation for the 
slow and fast magnetoacoustic modes has been found assuming 
evanescent-like perturbations in the coronal medium. Finally, the 
damping time of the fundamental oscillatory modes has been 
computed and the relevance of each non-adiabatic mechanism 
on the attenuation has been assessed. Next, we summarise the 
main conclusions of this work: 

1. Non-adiabatic effects are an efficient mechanism to obtain 
small ratios of the damping time to the period in the range of 
typically observed wavelengths of small-amplitude promi- 
nence oscillations. 

2. The mechanism responsible for the attenuation of oscilla- 
tions is different for each magnetoacoustic mode and de- 
pends on the wavenumber. 

3. The damping of the internal slow mode is dominated by 
prominence-related mechanisms, prominence radiation be- 
ing responsible for the attenuation in the observed wave- 
length range. Such as happens in the adiabatic case (see 
Paper I) the corona does not affect the slow mode at all, and 
these results are in perfect agreement with those for an iso- 
lated prominence slab. 

4. The attenuation of the fast mode in the observed wavelength 
range is governed by a combined effect of prominence radia- 
tion and coronal conduction. The presence of the corona is of 
paramount importance to explain the behaviour of the damp- 
ing time for small wavenumbers within the observed range of 
wavelengths. Non-adiabatic mechanisms in both the promi- 
nence and the corona are significant because the fast mode 
achieves large amplitudes in both regions. 

5. Since the external slow mode is principally supported by 
the corona, its damping time is entirely governed by coro- 
nal mechanisms, coronal conduction being the dominant one 
in the observed wavelength range. 

6. The consideration of different optical thicknesses for the 
prominence plasma causes an important variation of the 
damping time of the internal slow and fast modes in the ob- 
served wavelength range. Hence a precise knowledge of the 
radiative processes of prominence plasmas is needed to ob- 
tain more realistic theoretical results. 

7. The heating scenario has a negligible effect on the damp- 
ing time of all solutions in the observed wavelength range. 
Depending on the scenario considered, thermal instabilities 
can appear for small values of the wavenumber, in which 
coronal radiation dominates. 

8. The width of the prominence slab does not affect the results 
for both internal and external slow modes. However, fast 
modes are less attenuated in the range of observed wave- 
lengths when thinner slabs or filaments threads are consid- 
ered. 

Taking into account the results in the observed range of 
wavelengths, one can conclude that radiative effects of the 
prominence plasma are responsible for the attenuation of the in- 
ternal slow modes, which can be connected with intermediate- 
and long-period prominence oscillations, whereas a combined 
effect of prominence radiation and coronal thermal conduction 
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Fig. 13. Period (left), damping time (centre) and ratio of the damping time to the period (right) versus k x for the fundamental 
kink oscillatory modes: internal slow (top panels), fast (mid panels) and external slow (bottom panels). Solid lines correspond 
to z p = 3000 km whereas dotted lines correspond to z p = 30 km. Prominence (1) radiation conditions have been taken for the 
prominence plasma and the heating scenario is given by a = b = 0. 



governs the damping of fast modes, whose periods are compati- 
ble with those of short-period oscillations. 
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iBo dv z 
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Bo** 
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Appendix A: Expressions for the perturbations 

Combining Eqs. d7]l-(fT2]>. one can obtain the expressions for the 
perturbed quantities as functions of v z and its derivative 



Vr = 



Pi = 



P\ 



-ik x A 2 dv z 
u 2 -k 2 h 2 dp 
iu>po dv z 
oj 2 - k 2 A 2 dz ' 

iajpoA 2 dv z 
oj 2 - k\K 2 dz ' 



(A.l) 
(A.2) 
(A.3) 



Now, we write the expressions for the perturbations to the mag- 
netic pressure, p\ m , and the total pressure, p\i, 



Bo _ 'PQ v a dv z 

Pirn — oix ;r~ 

yu a) az 



Pit = P\+ P\ n 



ip (w 2 - k 2 x v\ ) dVz 



uk 2 



dz 



(A.7) 
(A.8) 



In the limit A — » c s (i.e. in the absence of conduction, radiation 
losses and heating), all the expressions reduce to those corre- 
sponding to the adiabatic case. 
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Appendix B: Approximate dispersion relation for 
the internal slow modes 

Internal slow modes are almost non-dispersive and for adia- 
batic perturbations a good approximation for the frequency is 
to * c sp k x , c sp being the prominence sound speed. In the non- 
adiabatic case, we can consider the equivalence between c s and 
A to propose oj as A p k x as an approximate dispersion relation. 
Taking into account Eq. (15[ for A, the approximate dispersion 
relation for the internal slow modes is a third order polynomial 
in a>, 



w 3 - iSco 2 - k 2 c 2 co + i—ml = 0, 

y 

with 

3K - (y - 1) (k|| p ^ + _ w pp) . 
S - (y- 1) {k\\ v k\ + oj Tp ) , 



(B.l) 

(B.2) 
(B.3) 



K\\p - % 



Pp 



In Fig. IB. II a comparison between the exact and approximate 
solutions is displayed and a perfect agreement is seen. 
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for the sausage modes. Quantities and S are given by 
Eqs. (|B.2t and (|B.3l l, respectively. 
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Fig. B.l. Period (left) and damping time (right) versus the lon- 
gitudinal wavenumber for the fundamental internal slow kink 
mode. The solid line corresponds to the exact solution and the di- 
amonds correspond to the approximate solution. Prominence (1) 
parameters and a — b = have been used in the computations. 



Appendix C: Dispersion relation for an isolated slab 

We here deduce a polynomial dispersion relation for the magne- 
toacoustic normal modes of a slab with a longitudinal magnetic 
field. Taking Eqs. (l27T i and (l28T l as the dispersion relations for 
the kink and sausage modes, respectively, one can replace k z and 
A with their correspondent expressions (Eqs. lfT4l - fT31 ), and the 
following fifth order polynomial equation is found, 

u 5 - iSaj 4 - (v 2 A + c 2 ) <K 2 o? + i [vlS + rfj- J 7CV 

+ v^clkl^tj - Uftv 2 A c 2 k 2 'K 2 = 0, (C.l) 
with 

+ ( n + l/2)V (n = ,l,2,...), 



for the kink modes, and 

7C 2 = ^ + !^ („= 1,2,3,...), 



